Theory of single spin inelastic tunneling spectroscopy. 
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Recent work shows that inelastic electron scanning tunneling microscope (STM) probes the el- 
ementary spin excitations of a single and a few magnetic atoms in a thin insulating layer. Here 
I show that these new type of spectroscopy is described using a phenomenological spin-assisted 
tunneling Hamiltonian. Within this formalism, the inelastic dl jdV lineshape is related to the spin 
spectral weight of the probed magnetic atom. This accounts for the spin selection rules observed 
experimentally. The theory agrees well with existing STM experiments for single Fe and Mn atoms 
as well as linear chains a few Mn atoms. The magnetic anisotropy in the inelastic dl/dV and the 
marked odd-even N effects are accounted for by the theory. 
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Electron tunneling is one of the central themes in 
condensed matter physics. It lies behind fundamen- 
tal phenomena like, Josephson effect^ and tunneling 
magnetoresistance 2 -, and provides an extremely versa- 
tile spectroscopic tool, both in tunneling junction^ and 
STM geometries 5 . The use of inelastic electron tunnel- 
ing to determine the vibration spectra of ensembles of 
molecules inside tunnel barriers goes back to the semi- 
nal work of Jaklevik and Lambe^. They observed steps 
in the differential conductance dl/dV curve at particular 
values of the bias voltages which matched the vibrational 
energy spectra of different molecules. This led to the no- 
tion of inelastic assisted tunneling^: an electron could 
tunnel across the barrier giving away its excess energy 
eV to create an elementary excitation. In this frame- 
work, as eV increases, new inelastic transport channels 
open, resulting in steps in the dl/dV curve. 

With the advent of the STM, it has been possible to 
downscale the technique of inelastic tunneling vibrational 
spectroscopy to the single molecule leve 1 !, a possibility 
anticipated in the early days of STM&. In a series of 
striking experiments ^ 10 ' 11 Heinrich et al have used in- 
elastic STM spectroscopy to probe the spin flip excita- 
tions of a single and a few transition metal atoms in a 
surface by means of STM spectroscopy. They have mea- 
sured the single Mn atom Zeeman gap2iL, the collective 
spin excitations of chains of up to 10 Mn atomsi - and 
the spin flip transitions within the ground state mani- 
fold of a single iron atom, split due to the single atom 
magnetic anisotropy 11 . Analogously, Xue et al. to have 
used inelastic STM spectroscopy to probe the spin ex- 
citations of one and a few Cobalt Phthalocyanines and 
have measured their exchange coupling^. 

Therefore, spin assisted inelastic tunneling spec- 
troscopy (SITS) provides a direct measurement of 
the spin dynamics of one or a few magnetic atoms, 
and complements spin polarised scanning tunneling 
spectroscopy^ which is sensitive to the average rela- 
tive orientation of the magnetic moments of tip and sur- 
face {friT and firs respectively). SITS permits to mea- 
sure energy scales like interatomic exchange J 10 ' 12 , g 
factor— i 11 ' 12 and magnetic anisotropy tensor—. These 
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FIG. 1: (Color online), (a) Scheme of the experimental setup 
(b) Elastic tunneling process, (c) Inelastic spin assisted tun- 
neling process, (d) Inelastic spin flip tunneling process. 



quantities determine the spin dynamics of the mag- 
netic atoms(s). In contrast to the case of vibrational 
spectroscopy^, the physical origin of coupling between 
the transport electrons and the local spins which makes 
SITS possible is not clear—. Hirjibehendin et al men- 
tion two possibilities^, exchange or dipolar coupling, 
although the former encompasses a variety of different 
mechanisms, like direct, kinetic, etc. The coupling must 
account for a number of experimental observations. The 
height of the steps in the dl/dV scales like the sum of 
the squares of the matrix elements between the initial 
and final states of the operators S a with a = x,y, z, the 
spin of the atom probed by the SITS 11 . This is related 
to the selection rule for the change of the spin of the 
local spin SS Z = ±1, 10 i 12 . Therefore, there is a rela- 
tion between the inelastic current and the spin spectral 
weight S aa (u>) of the magnetic atom(s). In this paper I 
show that an effective spin-assisted tunneling Hamilto- 
nian naturally explains the relation between the inelastic 
current and the spin spectral weight <S aa (o>) and accounts 
for the main experimental findings. 
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The experimental system consists of an insulating thin 
layer deposited on a metallic surface (see fig. la). Mag- 
netic atoms lie in the insulating layer and are probed 
by a STM . A natural model for this system would thus 
feature 3 types of fermion operators, tip, surface and in- 
sulating layer, plus the spin operator of the magnetic 
atoms. In such an approach, the current, evaluated to 
lowest order in the tunneling coupling, is related to to 
the spectral function of the transport electron in the in- 
sulating layer interacting with the local spins^. This is 
different from the experimental findings described above. 
In particular, the conductance so evaluated would have 
a Coulomb Blockade gap unless the central region is in a 
mixed valence point, but there dl/dV curve has peaks 
and not steps, and their location depends on the ex- 
change coupling between the local spin and the transport 
electrons^. 

Thus, within the 3 fermion approach the low bias steps 
found experimentally must arise from higher order cotun- 
neling processes^. Here I adopt a simpler approach using 
a phenomenological Hamiltonian with two types of elec- 
trons (tip and surface) , with an effective spin flip assisted 
tunneling termii: 

H = Hti p + H S ui +Hs + Htun (1) 

The first three terms describe the electrons in the tip 
(Wtip = J2k,a £ka f kcr a krT ) and in the surface (H SUT = 
a e kbp,jbp a ). The Hamiltonian of the central spin(s) 
is the sum of the intraatomic Hos(i) terms and the spin- 
spin couplings: : 

H S = £>os« + \ J2jab(i,j)S a (i) ■ S b (j) (2) 

The eigen-energies and eigenstates of Hs are denoted by 
= Em and \M). The tunneling terms can be written as: 

Htun= E T «( tt ')yS a (l)(4 5 W+h.c.) (3) 

kk'aa'a 

where t q and S a (l) are the Pauli matrices and the spin 
operators of the spin (1) in the central region for a — 
a = x,y, z, and the unit matrix for a = 0. Because of the 
short-range nature of exchange interaction and tunneling 
processes from the tip, we assume that only one spin 
S a (l) is assisting the tunneling. This term describes the 
tunneling of electrons between tip and surface assisted 
by exchange interaction with the spin in the middle. A 
similar term has been used by different authors^. 

Eq. ([3]) describes both spin assisted (a = x, y, z) and 
conventional tunneling a — 0. To lowest order in TLtun, 
the current has three conributions: (i) a central-spin in- 
dependent Tq, (ii) a crossed contributions proportional 
to ToT a mT,s • {iS'(l)) , and (iii) a spin flip contribution T 2 
that features the spin spectral weight and is responsible of 
the steps in the dl/dV curves. In the case of magnetized 
tip and sample the Tq contribution depends on mj- • mg 



which makes possible the SP-STS spectroscop y 1 ^ 19 . In 
this work rag = mj = 0. The spin-flip contribution 
arises from the balance between electrons tunneling from 
tip to surface and backii 

/ = E w) ( E «m s - E «,m t ] (4) 

M \ p,ff k,a I 

where P(M) is the equilibrium occupation of the M 
state, rip ,s is the occupation function of the tip and sur- 
face and 7 are the tunneling rates associated to the spin- 
flip assisted tunneling Hamiltonian: 

7 p T m 5 = E \T a (pp')\ 2 \(M\S a (l)\M')\ 2 x 

p',M' ,a 

X (l - rip,) S (e p > + e M > - e p - e M ) (5) 

This expression gives the lifetime of a product state 
with an electron in the state p of the tip and the mag- 
netic atom(s) in state M due to a spin flip assisted 
tunneling of the electron to the surface. Importantly, 
this equation relates current to the spin matrix ele- 
ments, \(M\S a \M')\ 2 ee \S MiM ,\ 2 , as reported in the 
experiment o 10 ^ 11 ' 12 . 

If the coupling between transport electrons and spins is 
rotationally invariant, T a [k, k') = T$ is the same for a 

y, z. The dependence of Tg on the momentum indexes 
can be neglected 2 ^. We take n T (e) = /(e) and n s (e) = 
/(e + eV), where / is the Fermi function. The sum over 
momenta leads to an integral over energies featuring the 
density of states of tip and surface, pt(e) and ps(e) which 
are assumed to be flat in the neighbourhood of the Fermi 
energy. We arrive to 

/oo />oo 
/ S aa (e-e')F(e 1 e\eV)dedt' (Q) 

where eGs = 1sPt(^f)ps{£f) and 
F(e, e',u) = [/(e) (1 - f(e' + w)) - /(e + W ) (1 - f(e'))} 
and 

S aa (Lu) ee Pm\SI UM 'W\ 2 5 (w + e M - e M >) (7) 

M.M' 

is the spin spectral weight. 

From the formal point of view, equation ((6]) is one of 
the main results of this paper. It relates the inelastic 
current to the spin spectral weight of the magnetic atom 
probed by the STM. It shows that two types of spin as- 
sisted tunneling processes contribute to the current. If we 
choose z as the quantization axis, the a = x, y terms in- 
volve spin exchange between the transport electron and 
the magnetic atom. These are the spin flip terms. In 
contrast, the a — z term conserves the spin of both car- 
rier and atom. This dichotomy is absent in the case of 
vibrational inelastic spectroscopy. 
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One of the integrals in eq. ([6]) is done using 
the delta functions in eq. ([7]) and the other using 

Ho t/( £ ) (1 " /( £ + de = T^F^' Thc total inelas " 
tic current is thus written as 



I = 



M,M'a 



PM\S a M<M ,(l)\ 2 i(eV,A M ,, M ) 



with A 



M'M 



E 



M' 



i(eV,A) = G s 



Em and 



eV- A 



eF + A 



j _ e -/3( e y-A) 



1 



(8) 



(9) 



is the current associated to a single inelastic channel with 
energy A. Equations © and © are the magnetic ana- 
log of the vibrational inelastic tunneling spectroscopy 6 
in which the dipole spectral weight of the molecular vi- 
brations is replaced by the spin spectral weight of the 
magnetic atoms 

Now the validity of eqs. (|6I8[) and the spin assisted 
tunneling term Q is verified by comparing their predic- 
tions with the experimental results. The case of tun- 
neling through a single Fe atom in CuN/Cu system^ is 
considered first. Following that reference, the standard 
single spin Hamiltonian reads 



DSl + E{S 2 X - Si) + g^ B B ■ S 



(10) 



with S = 2 adequate for Fe 2+ and D = —1.55 meV and 
E = 0.35 meV 11 . Since D >> E approximate analyti- 
cal expressions yield the ground state energy (for B = 0) 
E = -AD and the excitations &M'fl , 3D-6E, 3D + 
6E, 4D (see fig. 2c) . The ground (first excited) state is 
made mainly (only) with M z — ±2 . The dl/dV curves 
obtained from eq. ([5]) and the exact solution of Hamil- 
tonian (|10p for different intensities and orientations of 
the applied magnetic field, evaluated for fc^T = 0.5K, 
are shown in fig.2(a,b). The a (b) panel corresponds to 
field parallel to z (x). At zero field the dl/dV curves 
shows three steps corresponding to the excitations to the 
first, second and third excited states. The transition to 
the fourth state is forbidden ( Y, a =x,y,z l(°l 5 a|4)| 2 = 0). 
The prominent — > 1 transition comes from the spin- 
conserving channel a = z, where as the — ► 2 and 0^3 
transitions come from the spin flip channels a = y and 
a = x , respectively. Interestingly, the energy differ- 
ence between these two transitions is exactly equal to 6E. 
Thus, this parameter can be red from the experimental 
data. Their evolution as a function of the intensity and 
orientation of the magnetic field give good account of the 
main observed experimental features 11 . In particular, the 
conductance shows a significant magnetic anisotropy, re- 
lated to that of the iron atom in this surface. 

The theory also accounts for more complicated exper- 
imental situations where electrons tunnel through one 
magnetic atom which is exchanged coupled to others. 
This is the case of linear chains of N Mn atoms deposited 
on a CuN/Cu surface, N going from 1 to 10. Mn+ 2 has 
5 = 5/2 and very weak magnetic anisotropy. The Mn-Mn 



14 
12 

„ 10 

3. 8 

> 





- 5T 






- 3T 






-1T~l 






otI 




x _ 




I 





14 
12 
„ 10 

A 8 - 



4 h 
2 



\ w L 




1 


fr 











-5 5 
V (meV) 



-5 5 

V (meV) 



FIG. 2: (Color online), (a) Differential conductance for single 
Fe atom, with magnetic field along the z axis.(b) The same 
with magnetic field along the x axis, (c) Scheme of the 5 = 
energy levels of eq. (|10p and the tunneling induced transitions. 



coupling is approximated by a spin-rotational invariant 
first-neighbour Heisenberg coupling. The model reads: 

n S =J2 D iSltt) + WBB ■ S{i) + jJ2Sa{*) ■ S a {l + 1) 



where the sum in the last term runs from i = 1 to N — 1. 
Results are shown up to N — 4 spins for which the 
Hilbert space has 6 4 states. The experiments^ can be 
modelled taking J ~ 6 meV, much larger than the single 
Mn D. Thus, D is a weak perturbation of the Heisen- 
berg model, whose eigenstates can be labelled with the 
total spin S and its third component S z . Since j > 
the coupling is antiferromagnetic. Thus, even N chains 
have S — ground states whereas odd N chains have 
degenerate ground states with multiplicity 6 weakly split 
by the anisotropy term D. Thus, the lowest energy step 
occurs is related to D in odd N chains and to J in even 
N chains. 

The calculations (see figure ([3])) account both for the 
difference between odd and even N chains and the in- 
crease of the anisotropy related step of odd-iV chains ex- 
perimentally observed dl/dV^S-.The first excitation for 
N = 1 and N = 3 correspond to the transitions from the 
ground state doublet S z = ±5/2 to the first excited state 
with S z = ±3/2 through the spin-flip channels a = x,y. 
Application of a magnetic field increases the splitting be- 
tween these states and shifts the step towards higher en- 
ergy, as seen in the experiment^ and well captured by 
the model (fig 3b). 

The first excited states of the dimer have 3-fold degen- 
eracy weakly split into a low energy doublet and a higher 
energy singlet because of the single atom anisotropy. 
This fine structure splitting is 6AD. The ground state 
singlet is connected to the excited state singlet via a spin- 
conserving transition and to the excited state doublet via 
spin flip transition. Figures 3c and 3d show how this 
fine structure evolves as a magnetic field is applied. The 
model accounts for the deviation from the Zeeman split- 
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FIG. 3: (Color online), (a) dl/dV for 5 = 0, k B T = 0.6K 
for chains of N Mn atoms, N = 1,2,3,4. (b) Evolution the 
N — 1 low energy step as a function of B. (c) Evolution of 
the first step for N = 2 as a function of B. (d) Evolution of 
the corresponding exciations. 

ting observed experimentally 10 , with D — 0.08 meV and 
J = 5.89meV. Interestingly, the sign and magnitude of 
D for the dimer is different than that of the monomer. 

The numerical simulation also accounts for the red 
shift of the first step for even- TV chains as N increases^. 
This result can be rationalized using the Lieb, Shulz, 
Mattis theoremSi which states that the energy difference 
of the lowest energy excitation of linear chains with half- 
integer spin and the singlet ground state (even TV) must 
is bounded by 27r ^ s . Thus, the excitation energy is 
expected to decay as 1/N. 



In chains with N — 3 or more atoms the amplitude of 
the spin excitations, and the height of the steps in the 
dl/dV curves thereby, can vary from atom to atom. In 
figure 3a the dl/dV curve is shown for the spin-assisted 
current recorded on top of the central atom and one of 
the side atoms. As reported in reference fllll ) the location 
of the steps is the same. However, the intensity is clearly 
different. In the spectrum recorded in the side atom two 
steps are seen, corresponding to transitions to the first 
5 = 3/2 and second 5 = 7/2 excited states whereas the 
spectrum recorded in the central atom only the latter is 
seen. Thus, the spin-assisted tunneling spectroscopy can 
be used to map the amplitude of the spin excitations. 

In summary, the experiments of single spin inelastic 
tunneling spectroscopy imply the existence of a spin- 
assisted tunneling mechanism that couples electrons in 
the tip and the surface to the local spin ([3|). This term 
naturally leads to an expression for the current that in- 
volves the spin spectral weight S aa (uj) of the magnetic 
atom that assist the tunnel process and to the spin se- 
lection rules observed experimentally 10 ' 12 . Both the mi- 
croscopic origin of the spin assisted tunneling mecha- 
nism, and its connection to atomic scale current driven 
spin torque deserve additional work. The spin-assisted 
tunneling Hamiltonian is a genuine many-body process, 
since it involves four fermions, unlike ordinary tunneling. 
It could arise from the exchange part of the electron- 
electron repulsion in the Hamiltonian— and it could be 
a kinetic exchange term 2 ^. 

I acknowledge fruitful discussions with R. Aguado, J. 
J. Palacios, F. Delgado, C. Untiedt and C. Hirjibehedin. 
This work has been financially supported by MEC-Spain 
(Grant MAT07 ) and by Consolider CSD2007-0010. 

Note added:. During the completion of this work, a 
preprint with related results has been posted 2 ^. 
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